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INTRODUCTION 


The development and application of composite materials to aerospace vehicle 
structures which began in the mid 1960’s has now progressed to the point where what 
can be considered entire airframes are being designed and built using composites. At 
least two systems intended for production, the McDonnell -Douglas AV-8B and the 
Bell -Boeing V-22 aircraft are cases in point. Significantly, both are VTOL systems, 
in which empty -weight is an especially potent design variable. At the same time, 
certain aircraft and spacecraft components are being routinely designed and built using 
composites because either performance or economics require such use. At the low end 
of the temperature scale , tail surfaces provide one example of such applications 
(aircraft centers of gravity historically tend to fall too far aft); at the high end, 
reentry body lifting surface leading edges provide another example. 

As such applications have become more wide-spread, the initial developments and 
use of composites based on thermoset resin matrices are increasingly making way for 
those incorporating thermoplastics - for a number of reasons . Probably foremost 
among these are higher damage tolerance and wider allowable operating temperature 
ranges. Still more advanced concepts for such projects as the "Orient Express", 
National AeroSpacePlane (nasp), and the so called Advanced Tactical Fighter and 
Attack aircraft call for higher temperature applications of composites than have ever 
been sought before for airframes . in fact, missions which call for hypersonic fighter 
aircraft and missies, ultra -manueverable and/or VTOL fighters which make use of thrust 
vectoring nozzles, and transatmospheric vehicles which routinely leave the atmosphere, 
orbit, reenter, land and are capable of rapid turn-around for repeated reuse, are 
blurring some of the traditional distinctions between engine materials and structures 
and airframe materials and structures . For many of these concepts , higher 
temperature capable composite structural materials - over ranges which are rather 
continuous from ~ 400° to ~ 3000°f - must be considered, not just desireable, but 
enabling technologies. 

NASA and AFOSR continue to play leading roles in establishing the technology 
base required to realize the full promise of composites in sophisticated aerospace 
structures. This is being done through support of programs of fundamental research 
into composite materials and structures and the means by which they can be 
successfully applied in design and manufacture . RPI's program has been funded by 
NASA and AFOSR as part of an over -all university program in composites. Over the 



3 


roughly ten years in which it has been underway , its purpose has been to develop 
critical advanced technology in the areas of physical properties , structural concepts 
and analysis, manufacturing, reliability and life prediction. Specific goals have 
changed as the state of composite materials and structures art has developed. In the 
early years strictly low temperature airframe applications were of interest, and major 
efforts were expended in establishing new structural design concepts and in exploring 
low cost, innovative fabrication techniques. More recently, such research has given 
way to the pressing need to deal with the problems of higher operating temperatures. 
Furthermore, in this era, only the most fundamental aspects seem appropriate. 

The overall concept of RPI's program continues to be unusual for a university in 
several important aspects. First, the nature of the program has been comprehensive. 
We continue to probe to great depth in a relatively few, well chosen areas of 
investigation which, taken together, provide coverage of a wide spectrum of composite 
materials and structures issues. Although we have dropped many projects investigating 
the behavior of generic structural elements, fabrication science and technology and 
applicable, generally useful computer methodology developments, the expansion of 
temperature ranges has added a new dimension to the spectrum of composite issues. 
It is clear that as new directions are added to RPI's program, particularly with 
limited funds and personnel, older ones must be dropped. In this period, renewed 
emphasis was placed on the more fundamental issues associated with relatively 
little -explored areas of resin matrix composites and with the newer constituent 
materials: those fibers and matrices capable of higher temperatures. A number of 

studies involving the older resin matrix manufacturing processes , directionally 
solidified eutectics, and edge -initiated delamination failures have all been phased out 
over the past two years. 

Second, interactions among faculty contributing to program objectives is on a day 
to day basis without regard to organizational lines. These contributors are a group 
wider than that supported under the project. Program management is largely at the 
working level, and administrative, scientific and technical decisions are made, for the 
most part , independent of considerations normally associated with academic 
departments . This kind of involvement includes faculty , staff and students from 
chemistry, civil engineering, materials engineering, aeronautical engineering, 
mechanical engineering, and mechanics, depending on the flow of the research. 

Both of these characteristics of the NASA/AFOSR program of research in 
composite materials and structures foster the kinds of fundamental advances whch are 
triggered by insights into aspects beyond the narrow confines of an individual 
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discipline. This is often sought in many fields at a university, but seldom achieved. 

A third aspect is increasing the interaction between appropriate members of 
NASA's staff of Research Center scientists and engineers and those active in the 
program at RPI . This has required identification of individual researchers within 
NASA centers whose areas of interest, specialization and active investigation are in 
some way related to those of RPI faculty supported under the subject grant . A 
program of active interchange is then encouraged and the means by which such 
interaction can be fostered is sought. Benefits which result from this increased 

communication include a clearer window to directions in academia for NASA 
researchers; opportunities to profit from NASA experience, expertise and facilities for 
the faculty and students so involved; and an additional channel for cross -fertilization 
across NASA Research Center missions through the campus program. Finally, 

collaboration among RPI investigators is encouraged through management mechanisms; 
for example, asking faculty whose research promises to be synergistic to propose to 
the program’s Budget Advisory Committee jointly. 

In short, the NASA/AFOSR Composites Aircraft Program is a multi-faceted 
program planned and managed so that scientists and engineers in a number of pertinent 
disciplines at RPI will interact, both among themselves and with counterpart NASA 
Center researchers, to achieve its goals. Research in basic composition, 
characteristics and processing science of composite materials and their constituents 
has been planned each year, with the guidance of NASA and AFOSR technical monitors 
and Research Center engineers and scientists , to address the most pressing and 
promising aspects of composites in that particular era. In the current period, for 

example, issues related to the fabrication of non -resin matrix composites and the 

micro, mezzo and macromechanics of thermoplastic and metal matrix composites have 
been emphasized. 

In the following sections, more detailed descriptions of the progress achieved in 
the various component parts of this comprehensive program are presented. 


PART II 


RESEARCH 
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RESEARCH 


A. THE EFFECTS OF CHEMICAL VAPOR DEPOSITION AND THERMAL TREATMENTS ON 
THE PROPERTIES OF PITCH -BASED CARBON FIBER 

Sr. Investigator: R. J. Diefendorf 


INTRODUCTION 

chemical vapor deposition is one of the earliest means of processing composite 
materials and their constituents. It is of interest for carbon because of its potential 
for both increasing fiber properties and as a means of forming carbon/carbon 
composites . The purpose of this research is to determine the effects of carbon 
layers , established by chemical vapor deposition ( CVD ) , on the mechanical properties 
of pitch -based carbon fiber. One specific question is the extent to which a carbon 
coating can fill or "heal” those surface flaws which contribute to low stress failures. 

STATUS 

Work performed in the last reporting period indicated that moduli and possibly 
strength, as well as interfacial bonding could be altered with the application of a 
carbon coating. Another area of study was the effect on coating structure of different 
precursor gases. Initial work was also begun to determine if, by varying deposition 
parameters, fiber structure could be altered to improve its overall performance. The 
intent, of course, is to incorporate these findings, if successful, into carbon/carbon 
production techniques to produce a much higher performance material than is presently 
possible . 

PROGRESS DURING THE REPORTING PERIOD 

Furnace parameters were varied in deposition experiments in attempts to achieve 
the desired penetrative ability of the deposition species . Lower temperatures and 
pressures and shorter gas contact times yielded the best penetration results . The 
ranges of parameters used are listed in Table I-A-l, below: 
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T able ll-A-1 
Furnace Parameters 

Temperature 800-1600C 

Pressure 0.2-200 Torr 

Flow 50-6500 ccn/ e 

Time 0.16-10 hours 

Gases: Methane, Hydrogen, Naphthalene, and 
Dicyclopentadiene 

Pitch -based carbon fiber, in the form of a single filament or as a tow, was held in a 
graphite jig and used as a substrate for deposition. Tensile testing was performed to 
determine the mechanical properties SEM, X-Ray, and optical microscopy were used to 
analyze fracture surfaces, fiber surface topology, and deposition microstructure . 

A "sheath" effect occurred for the treatments in which a substantial coating was 
achieved. Fibers with apparently brittle and adherent coatings had reduced strength 
values. Coatings with little adherence and those that did not behave in a bittle 
manner either caused no adverse effects on properties or' improved fiber strength. 
Although coating thicknesses obtained during experimentation ranged from 
O to 2 jim, in many cases they were not discernible under SEM. In general, the fiber 
properties are not dependent upon the thickness of the deposited layer, but rather the 
parameters used to obtain the layer , i . e . , the structure of the layer . Lower carbon 
to hydrogen ratios and short gas contact times and/or high gas flow rates yielded in 
most cases, reduced strength and modulus values. 

Micros tructural changes resulting in improved fiber modulus were achieved when 
using dicyclopentadiene as a source gas at elevated temperatures of 160OC and in a 
3 -hour methane treatment at 1000C. More uniform preferred orientation properties 
appear to be a reason for the increased modulus values; however, the mechanisms for 
this improvement are still being studied. 

In other experiments to determine if the increased mobility of carbon atoms in the 
presence of hydrogen could lead to improved preferred orientation of the fiber 
structure, it was determined that a hydrogen environment, with or without the addition 
of an applied stress, has little effect on the fiber's structure or its mechanical 
properties at the temperatures investigated. 



8 


PLANS FOR THE UPCOMING PERIOD 

Studies to determine the actual mechanism responsible for the increased modulus 
values observed in certain treatments will be continued. Modulus can be altered by 
improving preferred orientation, by decreasing the interlayer compliance of the fiber, 
or by the addition of a highly oriented , adherent coating . The experiments will 
center on these effects. More work at the higher temperatures will be performed to 
observe the reproducibility of the "sheath" effect and its effect on strength of the 
fiber. The possibility of increased strengths with higher temperature CVD seems to be 
indicated . 

Finally, the effect of CVD coatings from different precursor gases also will be 
investigated to determine their effect on interfacial bond strength. By controlling the 
interface with a CVD coating, vast improvements in the properties of epoxy matrix and 
ceramic matrix composites could be achieved. 

PRESENTATIONS AND PUBLICATIONS BY PROF. R. J. DIEFENDORF ON THIS SUBJECT. 


"The Relationship of Structure to Properties in Carbon Fibers", presented at the Third 
Japan-U.S. Conference on Composite Materials, Tokyo, Japan, June 23-25, 1986. 

"The Chemical Vapor Deposition of Carbon Capillary Tubes II", presented at the 4th 
International Carbon Conference, Baden-Baden, Deutschen Keramischen Gesellschaft , 
West Germany , June 30 -July 4, 1986. 

"Carbon Fibers From Mesophase Pitch", presented at the Rockwell Science Center, 
Thousand Oaks, Ca. July 11, 1986. 

"Carbon Fibers", presented at SUNY- Buffalo, N.Y., April 2, 1987. 

"Chemical Vapor Deposition" , presented at Sohio , Niagara Falls , N . Y. , April 6 , 
1987. 

"Composite Processing", presented as the Sach's Memorial Lecture, Syracuse 
University, Syracuse, N.Y., April 21, 1987. 
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B. INELASTIC DEFORMATION OF METAL MATRIX LAMINATES 


Sr. Investigators E. Krempl 


INTRODUCTION 

Classical laminate theory, assuming uniformly or linearly distributed strains (for 
in-plane loading and bending, respectively), has served as a useful tool in the design 
of composites Cl, 21. This theory is presently limited to linear elastic behavior and 
cannot be employed when significant time dependence or plasticity occurs. Both such 
effects cam be encountered in metal matrix composite applications , especially when 
elevated temperature service is contemplated. 

STATUS 

During the past decade, the principal investigator and his students have examined 
the room temperature and elevated temperature behavior of engineering alloys using 
servocontrolled testing machines . The behavior of the alloys was found to be 
viscoplastic , even at ambient temperatures C3 through 71. Consequently, the theory 
of viscoplasticity based on overstress ( VBO) was developed in uniaxial C81 and in 
isotropic form C91. A full invariant, orthotropic version of this theory has been 
established using tensor function representation theorems C101. A simplified version 
of this theory was derived in [ill. This simplified version retains the essential 
features of the original theory including the absence of a yield surface associated with 
loading and unloading conditions and the existence of asymptotic solutions for constant 
strain (stress) rate loading. In the theory, the number of material functions is 

reduced as much as possible to facilitate applications. 

PROGRESS DURING REPORTING PERIOD 

Lamina Behavior 

The simplified theory of [ill was specialized for the plane stress case so as to 
represent the state of stress in a single ply. Details cam be found in a forthcoming 
report 1121. 

To demonstrate the usefulness of the theory, it was applied to test data reported 
by Kreider and Prewo [131 and reproduced here as Figure II-B-1. In this figure, 
uniaxial stress -strain diagrams in various directions are shown for a single ply of 
BORSIC( 100-jim)/AL-6061-T6 metal matrix composite. Since no test results for the 
viscous properties of this metal matrix composite were provided in C131, the viscous 





Figure II-B-1 

train Curves for 6061 Aluminum Reinforced with 100-ym Borsic Fibers 
Tested at Indicated Angles to Fiber Axis, From Kreider and Prewo [13]. 
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properties of AL-6061-T6 alloy given in C91 were used. The matrix in C131 and the 
alloy in C91 have the same designation and heat treatment, but it is recognized that 
in situ and neat metal properties can be different. 

Figure II- B -2 shows the correlation provided by the orthotropic VBO theory using 
a particular set of material constants. The curves represent the numerical integration 
of the system of differential equations under constant strain rate. 

All integrations were performed on an IBM AT personal computer using the routine 
DGEAR for the solution of stiff nonlinear differential equations. 

Once the material functions and constants are determined, the set of differential 
equations can be used for other deformation histories. As an example, the lamina 
behaviors predicted by the theory under cyclic loading are shown in Figures II -B -3 
and 4 . The loops shown here are closed after one cycle , demonstrating the 
consequences of a theory incorporating cyclic neutral behavior. Figures Il-B-3a and 
3b show hysteresis loops with short time relaxation periods in directions relative to 
the fiber axis of 5° (Fig. II-B-3a) and 45° (Fig. II-B-3b), respectively. It can 
be seen that relaxation behavior is much more pronounced on loading than on 
unloading, and that it is more pronounced in Figure II-B-3b than 3a, due to the 
increased influence of the matrix (which is viscoplastic). Predictions of short-term 
creep and hysteresis are shown in Figures II-B-4a and 4b. Creep is also more 
pronounced during loading than unloading . No creep is observed during unloading 
where the material behavior is nearly elastic . 

No cyclic experiments were reported in C133, so that no cyclic comparisons are 
possible. The predicted hysteresis loops show the essential features of cyclic neutral 
behavior, however, and the time dependence is also sub expected based on isotropic 
material behavior at room temperatures C51. Until cyclic experiments with metal 
matrix composites become available, the features shown in Figures II -B- 3 and 4 will 
remain uncorroborated theoretical predictions. 

Laminate Behavior 

A laminate theory for metal matrix composites was then developed following the 
assumptions and methods of the classical laminate theory Cl, 21, and incorporating the 
orthotropic VBO theory. This will also be reported in C121. The equations were 
programmed for routine DGEAR integration of the stiff differential equations of the 
viscoplasticity theory based on overstress which result for cases representing the 
in -plane loading of symmetric but otherwise arbitrary lay-ups. 

Examples of the cyclic behavior predicted by this theory are given in Figures 
II -B -5 and 6. In each case, one cycle is computed with relaxation periods as in the 
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Figure II-B-3a 

Cyclic Loading of a 5° Ply at a Strain Range of 1%. Predicted by the VBO 
at a Strain Rate of 10~ s 1/s. At Points A, A', B, B' a 1000 s Relaxation 
Hold-Time is Introduced. Note that the Stress Relaxes Only at Points A and 
A' During Loading but no Relaxation is Observed at the Same Strain During 
Unloading, Points B and B'. 
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20. L 
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STRAIN 


Figure II-B-3b 

Cyclic Loading of a 5° Ply at a Strain Range of 1% in the 45° Direction. 
Due to Matrix Dominance, the Hyseresis Loop in the Same Strain Range is 
Much More Pronounced than in Fig. II-B-3a. Again the Relaxation Drop is 
Much More Pronounced at Points A, A' than at Points B, B'. 
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Figure II-B-4a 

Same as Fig. II-B-1 Except that a 200 s Creep Hold-Time is Introduced at 
Points A, A', B, and B'. Although the Magnitude of the Stress is Equal at 
Points A, A* and B, B', no Creep is Noticeable on this Graph at Points B, B 
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Figure II-B-5 r 

Behavior of a [±30]s Laminate Predicted Using VBO. A 1000 s Relaxation 
Hold-Time is Introduced. at Points A, A', B, and B'. The Difference of 
the Relaxation Behavior at Points A, A', B, and B' is Preserved. Strain 
Rate 10 -s 1/s. 
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-.5 


STRAIN 


Figure II-B-6 

Same as Fig. II-B-5 Except that the Laminate is a [0/90] s Lay-Up. Th 
Fiber Dominance in the Deformation is Evident by the Small Hysteresis 
Loop and the Insignificant Stress Drop During Relaxation. 
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lamina case. It is evident that the behavior of the [0/903 3 laminate in Figure 
II-B-6 , is quite different from the matrix dominated behavior of the [±303 s laminate 
shown in Figure n-B-5. 

As in the case of the cyclic ply computation, no experimental data are available 
for correlation with these laminates. These figures are, therefore, mostly intended to 
demonstrate the versatility of the proposed approach. Details will be reported in 
[123. 


PLANS FOR THE UPCOMING PERIOD 

The theory will be correlated with the results of other experiments available in 
the literature. A simple laminate theory for bending will be formulated and numerical 
experiments performed. 

PRESENTATIONS AND PUBLICATIONS BY PROF. E. KREMPL ON THIS SUBJECT 

"The Viscoplasticity Theory Based on Overstress Applied to the Modeling of a Nickel 
Base Super Alloy at 815*C", with Lu, H. a Yao, D., presented at the 3rd Symposium 
on Nonlinear Constitutive Relations for High Temperature Applications, Akron, Ohio, 
June 11-13, 1986. To appear in the Proceedings. 

"Cyclic Uniaxial and Biaxaiai Hardening of Type 304 Stainless Steel Modeled by the 
Viscoplasticity Theory Based on Overstress", with Yao, D., presented at the 3rd 
Symposium on Nonlinear Constitutive Relations for High Temperature Applications. 
Akron, Ohio. June 11-13, 1986. To appear in the Proceedings. 

"A Simplified Orthotropic Formulation of the viscoplasticity Theory Based on 
Overstress", with Sutcu, M, presented at the 3rd Symposium on Nonlinear Constitutive 
Relations for High Temperature Applications , Akron, Ohio, June 11-13, 1986. To 
appear in the Proceedings. 

"Biaxial Fatigue and Deformation Behavior of Graphite/Epoxy Composites", presented 
seminar. Department of Mechanical Engineering, University of Delaware, Newark, 
Delaware, September 19, 1986. 

"Biaxial Fatigue and Deformation Behavior of Graphite/Epoxy Composites", presented 
at Southwest Mechanics Lecture Series, The University of Oklahoma, Norman, 
Oklahoma j November 10, 1986. 

"Isotropic and Orthotropic Formulations of the Viscoplastic Theory Based on 
Overstress", presented at International Conference on Constitutive Laws for Engineering 
Materials! Theory and Applications, Tucson, Arizona, January 7, 1987. 
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C. ANALYSIS OF FATIGUE DAMAGE IN FIBROUS MMC LAMINATES 

Sr. Investigator: G. Dvorak 


INTRODUCTION 

The mechanism of fatigue deunage in Metal Matrix Composite (MMC) laminates can 
be described as follows: cracks nucleate and grow in individual plies of the laminate 

as a result of cyclic plastic straining of the matrix . If the plastic straining is 
terminated , e . g . , by ply cracking which allows a crack accommodation strain to 
replace the cyclic plastic strain, then the matrix, or more precisely, each matrix 
segment between cracks, returns to an elastic state. If this happens in all plies, the 
composite "shakes down", damage accumulation stops, a saturation damage state is 
reached. Our objective is to find the crack density in each ply which corresponds to 
the shake -down state of a laminate under applied cyclic load. 

STATUS 

This work is part of a continuing investigation of the mechanisms of fatigue 
damage in metal matrix fibrous laminates. Earlier experimental work C14J showed that 
B-Al laminates under cyclic loading may experience extensive matrix cracking that, in 
a saturation state, may cause a substantial reduction (~-50%) of overall stiffness and 
strength. The same effect was later observed with SiC-AJ laminates. We had also 
found that the extent of fatigue damage depends on applied load amplitude and that no 
damage occurs in laminates cycled within the shake -down range, in which there is no 
cyclic plastic straining in the matrix. Results obtained in the present work suggest 
that the extent of fatigue damage at a given load amplitude is determined by the 
requirement that the composite laminate must reach a shake -down state through damage 
accumulation. 

Progress During the Reporting Period 

The analysis conducted during the reporting period is based on a combination of 
techniques derived from plasticity analysis of MMC plates and from damage analysis of 
elastic composite laminates C153, C161 , C17J. For each given increment of load 
amplitude, one finds an increment in ply crack density that assures that the average 
stresses in the ligaments that remain between the cracks do not violate the yield 
condition for the ply in question. 

The analysis was performed in strain space, because all plies of the laminated 
plate experience identical strain magnitudes under in -plane loads. Plastic response of 


21 


the plies earn then be described, in part, with the help of relaxation surfaces. The 
ply remains elastic if strained within its relaxation surface. If cracks are added to 
the ply, the overall stiffness of the ply is reduced and, therefore, the relaxation 
surface expands for strains which cause the cracks to open. An example is shown in 
Figure U-C-l. In a cyclically loaded laminate, cracks are incrementally added to all 
plies in which the yield condition would be violated under current overall load. 
Figure II-C-2 shows an example of relaxation surfaces of a (0/90) 2S laminate which 
has reached a certain damage state. 

The modeling procedure was applied to severed B-Al laminates which were tested 
in an earlier experimental program. Stiffness changes caused in saturation damage 
state at various load amplitude levels were calculated and compared with experiments. 
The theoretical results were found to be in good agreement with experimental data. 
This is shown in Figure Il-C-3. 

PLANS FOR THE UPCOMING PERIOD 

In the 1987-88 research program we hope to develop new models of 
time -dependent deformation of fibrous metal matrix Composites at elevated 
temperatures. The approach will be based on our previous work in plasticity of MMC, 
and it will incorporate the creep properties of the fiber and matrix. In particular, 
we expect that the fibrous composite may deform in several distinct modes, depending 
on the applied state of stress , and on the elastic properties of the phases . 
Micromechanical analysis of these modes will be performed, and predictions of overall 
instantaneous response will be made in terms of phase properties , microstructural 
geometry, and previous deformation history. 


PRESENTATIONS AND PUBLICATIONS BY PROF. G. DVORAK ON THIS SUBJECT 

"Damage Mechanics of Composite Materials" , colloquium, Northwestern University, May 
2, 1986. 

"Fracture Mechanics of Metal Matrix Composites” , presented at ONR Workshop on 
Failure Mechanics, University of Maryland, May 12-13, 1986. 

"Analysis of Fatigue Cracking of Fibrous Metal Matrix Laminates" , Symposium on 
Advanced Composite Materials, General Electric Co., May 14, 1986. 

"Damage in Metal Matrix Composites", presented at SDIO/ONR Composites Consortium 
Program Review, Woods Hole, MA. , June 2-3, 1986. 

"Thermal Expansion of Elastic -Plastic Composite Materials", presented at 10th U.S. 
National Congress of Applied Mechanics, ASME, Austin, Texes, June 16-20, 1986. 
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Figure II-C-1 

Relaxation Surfaces of the 90° Layer at Different Crack Densities 
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PERCENT OF INITIAL ELASTIC MODULUS 
AFTER 2x 1 0 # CYCLES 


Figure II-C-3 

Change in Elastic Modulus of a B-Al Plate Related to Applied Stress 
Range. Comparison of Theoretical Predictions with Experimental Data 
reported by Dvorak and Johnson (1980) and Johnson (1979). 
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"Thermomechanical Couplings in Solids" and "Thennomechanical Deformation and 
Coupling in Elastic -Plastic Composite Materials", invited lecturer, IUTAM Symposium, 
Paris, France, September 1-5, 1986. 

"Plasticity of Composite Materials”, colloquium, Texas A&M University, November 11, 
1986. 

"Plasticity of Composite Materials", colloquium, Rice University, November 12, 1986. 

"Fatigue Damage Analysis in Metal Matrix Laminates”, invited lecturer, AS ME Winter 
Annual Meeting, Anaheim, CA, December 7-12, 1986. 

"Recent Developments in Plasticity of Fiber Metal Matrix Composites" , Mechanics 
Colloquium, Yale University, January 21, 1987. 

"Bimodal Plasticity Theory of Composite Materials”, Symposium in Memory of Aris 
Phillips, Gainesville, Florida, January 28-30, 1987. 

"Recent Developments in Plasticiity of Composite Materials" , seminar , University of 
California at Berkeley, February 23, 1987. 

"Cracks Approaching Interfaces: The Image Crack Method", ONR Workshop on 

Composite Materials - Interface Science, Leesburg, Virginia, March 11, 1987. 

"Analysis of Metal Matrix Composites for Spacecraft Applications", presented at ONR 
SDI Review, University of Maryland, March 30, 1987. 

"A Bimodal Plasticity Theory of Composite Materials", seminar. Brown University, April 
6, 1987. 

"Recent Developments in Composites Plasticity" , seminar , Lawrence Livermore 
Laboratories, Livermore, CA, April 23, 1987. 
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D. DELAMINATION FRACTURE TOUGHNESS IN THERMOPLASTIC MATRIX COMPOSITES 


Sr. Investigator: S. S. Sternstein 


INTRODUCTION 

The delamination fracture toughness of high performance composite laminates is of 
considerable importance in applications where some degree of out -of -plane loading will 
be experienced or in cases where in -plane compression loads must be supported. In 
the latter case , delamination severely limits the compression loads which can be safely 
carried. Many investigators also believe that delamination toughness plays a major 
role in damage tolerance with respect to planar impacts. 

Thermoplastic matrix composites have the potential for significant improvement in 
delamination fracture toughness. However, the relationship between neat matrix and in 
situ matrix behavior is not understood well , owing in large measure to the finite 
strain , nonlinear viscoplastic behavior which thermoplastic matrices are known to 
exhibit . The objective of this research is to investigate the parameters which 
influence such behavior. 

STATUS 

Several studies of the micromechanics of delamination failure as related to matrix 
properties have been initiated which are aimed at elucidating the matrix -related failure 
modes of thermoplastic matrix composites. Specifically, the following studies are in 
progress: 

1) Delamination fracture toughness tests aimed at relating the Double Cantilever 
Beam ( DCB) delamination compression buckling fracture toughness to ply microstructure 
and local deformation (crack path) patterns. 

2) A companion series of numerical analyses, to examine the DCB test itself 
using finite element analysis (FEA) and a nonlinear, rate and stress state dependent 
constitutive equation specifically developed for the polycarbonate matrix material used 
in Part 1. This study is joint with Professor M. Shephard who is responsible for the 
FEA computations. (See Part II -E of this report.) 

3) New experiments, initiated to examine delamination fatigue crack initiation 
mechanisms in several thermoplastic and thermoset matrix composites. 

4) Companion numerical analyses, again using advanced constitutive relationships 
and FEA for the Mode II delamination fatigue problem of Part 3. 

5) A related study, currently in progress, on the compression strength of 
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thermoplastic matrix composites. 

PROGRESS DURING THE REPORTING PERIOD 

This progress report will cover only Part 1 of the program, as listed above. Our 
intent here is to develop a self-consistent data set applicable to the delamination 
fracture toughness of polycarbonate matrix composites, with and without fiber sizing, 
and the resultant fracture path as examined by reflected light microscopy . In 
addition , electron microscopy studies were performed on the sized and unsized samples 
to determine the extent of fiber -matrix adhesion. 

Experimental Results 

Unidirectional, 12 ply laminates were prepared using high quality polycarbonate 
matrix prepreg prepared by NASA-Gangley . Both unsized and epoxy sized AS4 carbon 
fibers were used. In some, but not all, samples an additional film of polycarbonate 
of either 50|xm or 250fim thickness was inserted at the midplane. In ail cases a 
starter crack was introduced by using a Kapton film at one edge between the central 
plies of the laminate. Double cantilever beams were cut from the master sheet using 
standard techniques. Each sample was fitted with aluminum end -blocks containing holes 
for pivot pins through which the loading was applied. The end -blocks were epoxied to 
the samples. Samples were deformed in the DCB mode using a crosshead rate of 1 
cm/min. An LVDT (linear variable differential transformer) was used to measure 
overall displacement at each of the two ends of the double cantilever beam. 

Microscopy samples were obtained from the broken DCB samples and either potted 
for reflected light microscopy or gold-plated for scanning electron microscopy. 

A typical load vs. deflection curve which exhibits permanent deflection is shown 
in Figure ll-D-1. Permanent deflection is accounted for in the data reduction 
scheme. Most samples demonstrated at least seven crack jumps (non -planar crack 
propagation). The inset graph shows the sample compliance vs. the edge -measured , 
crack jump length. Correction for the plastic deformation offset of the load -unload 
curves was made. The compiance vs. crack jump curve was fit using a second order 
Chebyshev polynomial which was differentiated to give dc/da. 'Hie standard DCB 

p2 

formula was then used to compute G lc = — (dc/da) where P is the breaking 
load, w is the sample width and dc/da is the change in sample compliance with 
crack jump length. The results are given in tabular form in Table Il-D-l and 2 for 
the various samples. The values of Gj c reported are the averages for each sample 
excluding the first crack jump. The standard deviation for the values of Gj c from all 
the crack jumps obtained in a given sample are also shown (as a) . For the 51jxm 
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avg . among 
samples 


avg. among 
samples 


Table Il-D-l 


Polycarbonate - AS 4 Unsized 
DOB Fracture Toughness 

No Central Film Slum Film 254um Film 


G Ic 

A 

G IC 

A 

G Ic 

A 

(J/m 2 

) 

(J/m' 

S ) 

(J/m 2 ) 


859 

95 

1179 

120 

1564 

238 

914 

234 

1319 

136 

1495 

168 

1105 

268 

1549 

270 

1234 

86 

1235 

194 

990 

198 

1352 

77 

825 

142 

1482 

208 

952 

100 

1090 

138 

1298 

222 

1470 

216 

886 

159 

1124 

165 

1516 

158 



856 

268 

1516 

117 



1026 

210 

1650 

184 



1047 

92 

1729 

152 



920 

129 





1283 

926 



987 

162 

1173 

218 

1447 

222 


Table II -D -2 

Polycarbonate - AS 4 Epoxy Size 
DCB Fracture Toughness 


No Central 

Film 

Slum central 

Film 

254 tun Film 

G Ic 

(J/m 2 ) 

A 

G Ic 

(J/m 2 ) 

A 

G lc 

(J/m 

A 

2 ) 

2217 

236 

2129 

60 

1814 

144 

2197 

289 

2565 

96 

1916 

205 

2415 

278 

2000 

190 

1753 

332 

2273 

183 

2081 

184 

2131 

598 

2338 

248 

2362 

203 

2160 

363 

2209 

250 

2114 

287 

1920 

233 

2151 

231 

2331 

171 

*1733 

167 

1875 

298 

2023 

176 

1769 

146 

2470 

302 

2065 

173 

*2162 

119 

1910 

128 

2131 

151 

2236 

213 

2170 

202 

1949 

222 



2139 

425 





2135 

222 





2188 

163 

2159 

185 

1959 

195 


*See Figure II-D-6 and Figure 7 for optical micrographs of these samples 
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film, unsized sample which showed a standard deviation of 926 J/m 2 (See Table 
H-D-l), removed of the second and third crack jumps reduces the standard deviation 
to 48 J/m 2 . Similarly, die standard deviation of 598 J/m 2 for the 254fun film, 
epoxy-sized sample (See Table H-D-2), was reduced to 147 J/m 2 by removal of the 
second crack jump. 

It is particularly noteworthy that there was not a regular pattern of increasing or 
decreasing Gxc vs. crack jump number in any given sample. Also, the standard 
deviation of Gr. samples is not substantially different from that calculated from the 
results obtained within a sample . Other investigators have observed a regular 
increase of Gj c with each crack jump in a given sample and have attributed this to 
fiber bridging . We see no evidence of such an effect and speculate that fiber 

wavyness and/or fiber wash during lamination may have resulted in trends not 
supported by our data. In any event, this again illustrates the lack of agreement 
among laboratories on specific properties of composites. 

The epoxy-sized samples (Table II-D-2) show no effect of film thickness on 
measured Gj c . However, the 254/im film in the unsized samples (Table II-D-l) 
appears to have given a statistically significant higher Gj c them the other unsized 
samples. The fracture toughness of the epoxy -sized samples is about double that of 
the unsized samples. This suggests that there is a favorable influence of increased 
interfacial adhesion between matrix and fibers. Such a conclusion is borne out by the 
SEM photomicrographs in Figures n-D-2 and 3. The presence of ribbons of 

debonded polycarbonate matrix are clearly evident in Figure II-D-3 (unsized) Whereas 
the adhesion of matrix (grey) to the fibers is evident in Figure H-D-2 (sized). 
Higher magnifications of the same failure surfaces sure shown in Figures II -D -4 and 5. 

Optical micrographs are shown in Figures H-D-6 and 7 for two samples cut from 
the same master sheet, containing a centred film of 254|un thickness and with 

epoxy -sized fibers. One sample tested at the low end of the Gxc values for this 

sample batch, whereas the other tested high for Gjc- In Table n-D-2 , these samples 
are indicated with an asterisk. It is apparent that the sample with a low Gj c (Figure 
n-D-6) contains a delamination crack which very infrequently crossed the central 
tough film. Conversely, the tough sample (Figure H-D-7) shows a crack path which 
has jumped across the tough film resulting in a greater area production rate per unit 
crack jump. 

We offer no explanation for the determining factor whcih controls the crack 
jumping (or the lack of such jumping), but conjecture that residual stresses (thermal, 
etc . ) and local variations in microstructure resulting in variations of stiffness and 
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Figure II-D-4 

Epoxy Sized Sample, 254ym Film. Matrix Failure, Limited Debonding. 

G t = 2235 J/m 2 . 

Ic 



Figure II-D-5 

Unsized Sample, No Films. Fibers Debonding. G = 1233 J/m^. 
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Figure II-D-6 


Epoxy Sized Sample, 254ym Film. Failure Occurs Within Ply. 


Ic 


1733 J/m 2 


serted Film 


:ting Agent 





Figure II-D-7 
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strength play a role. Me are attempting to model such effects using PEA. in any 
event, it is clear that crack jumping (nonplanar crack propagation) may play a 
significant role in the statistical spread of fracture toughness values among ostensibly 
identical samples. 

PLANS FOR THE UPCOMING PERIOD 

In the next reporting period we plan to place major emphasis on the constitutive - 
microstructural modeling of failure properties in thermoplastic matrix composites and to 
investigate both Mode 1 and Mode 2 delamination problems using PEA methods. 


PRESENTATIONS AND PUBLICATIONS BY PROF. S. STERNSTEIN ON THIS SUBJECT 

"Thermoplastic Matrix Composites", seminar , University of Connecticut, Storrs, CT, 
March 8, 1987. 

"Matrix Dominated Mechanical Properties of Composites”, seminar, TRW Research 
Center, Los Angeles, CA. , April 10, 1987. 

"Mechanical Properties of Composites”, lecture, Alcoa Conference on Micromechanics, 
Hilton Head, SC, June 8, 1987. 

"Thermoplastic Matrix Composites”, lecture, Gordon Conference on Composites, Santa 
Barbara, CA., January 1987. 

”A Micrographic Study of Bending Failure in Five Thermoplastic/Carbon Fiber 
Composite Laminates”, with S. w. Yurgartis, accepted for publication by Journal of 
Materials Science, April 1987. 
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E. NUMERICAL INVESTIGATION OF THE MICROMECHANICS OF COMPOSITE FRACTURE 

Sr. Investigator: M. S. Shephard 


INTRODUCTION 

The behavior of and failure mechanisms in composite materials can usually be 
better understood if analytic/numerical analyses are used appropriately in support of 
an experimental program. The goal of this work is to provide nonlinear finite element 
analysis capabilities which , coupled with am experimental program , will provide 
increased understamding of composite behavior. Initial efforts consist of two phases : 
first, modeling thermoplastic composites, concentrating on the development of nonlinear 
time -dependent constitutive relations; amd second, incorporating them into an existing 
nonlinear finite element program that ham the capability of adding new material models. 
Professor S . Sternstein is providing the technical expertise required to develop 
constitutive equations that can realistically model the measured material behavior . 
(See Section H-D, of this report.) 

STATUS 

The currently available finite element amalysis codes which could provide a 
framework on which to build analyses for the nonlinear analysis of composite materials, 
do not contain the forms of constitutive relations and mixing models necessary to 
properly analyze these materials in the nonlinear range. Therefore, our efforts have 
been concentrated on the development and implementation of such models in a nonlinear 
finite element code; namely, ABAQUS . As indicated below, the development and 
implementation of these procedures is nearing completion . In addition , numerical 
studies of specific experimental configurations were initiated. 

PROGRESS DURING REPORTING PERIOD 

Efforts during the reporting period concentrated on developing a time -dependent 
constitutive relationship for the matrix material and combining this mathematical model 
with the periodic hexagonal array model of fibers in this matrix. 

Matrix Constitutive Relation 

A nonlinear viscoelastic material model has now been developed. Emphasis was 
placed on capturing the essential mechanical characteristics of thermoplastics , 
including rate dependence, stress component interactions and transient behavior. The 
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Eyring stress -biased barrier model is generalized to three dimensions such that the 
hydrostatic pressure correctly modifies the deviatoric response. 

The rheological model chosen to represent the viscoelastic response of 
thermoplastic matrix material is based on the following two assumptions: 

1) The linear -elastic hydrostatic component of the hydrostatic stress produces 
an elastic dilatational response represented by 

CT kk = 3Be kk 

Where B = Bulk modulus 
cr = stress (total) 
e = strain (total) 


2) The non-linear rate -dependent deviatoric component of strain, e, can 
be related to stress by assuming that the material is represented 
by a Maxwell element in series with a Voigt element, as follows 

(in which ( ' ) ^ “ ): 


e - 



S_ 

n i 



(e . 



S_ dt) 
o n i 


CEq.ll 


Where n = solid state viscosity 
S = deviatoric stress 
G = shear modulus 

Subscripts 1 and 2 define properties before and after 
onset of nonlinearity, respectively, 
e = e - Vstre-I 
S = ct- Vatre.I 

tre indicates the trace of the matrix 
and I = unit matrix 


For the dashpot in the Voigt element, an activated non -Newtonian viscosity term is 
used (I. M. Ward C191), which obeys the Eyring hyperbolic -sine flow law. Ey ring's 

model assumes that applied stress shifts the potential energy barriers decreasing the 
material's resistance towards deformation. 

In the Eyring model: e = S dashpot . e = K sinh(a S dashpot ) 
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AH 

Where: K = e Q e Rr 

v 

a = — 

RT 

and e 0 = a constant pre-expontial factor 

v = activation volume for the molecular event 
aH = the activation energy (assumed constant by Eyring) 
R = ideal gas constant 
T = absolute temperature 


Sherby and Dorn (i.M.Ward C19 3) had found that increasing stress decreases 
the activation energy; ah •» ah - where /x = parameter which controls the 

decrease of the activation energy with increasing hydrostatic stress and 

CT m = i a xx 


Sternstein and Ho (I.M. Ward C193) found that hydrostatic stress causes 
the onset of nonlinearity. These two ideas were combined to yield 

e = K-e^inh(aS dashpot ) CEq. 21 


where (3 — 

RT 


Substituting CEq. 23 into CEq. 13 led to the relation 


*ii = + + Re* 30 " 1 • sinh 

2G X n 




Sij - 2G 2 ( eij - ^ - 


Sii dt, 


CEq. 33 


The resulting equation was successfully used by Cessna and Sternstein C203 
to characterize polymer deformation at the leading edge of a crack tip. 

Since a displacement -based Finite Element Method was anticipated, this 
equation was inverted to the form: 

S = f(e) 

Then, using deviatoric components instead of total stresses, the final form 
of the general, 3-D nonlinear viscoelastic constitutive equation was obtained; 
namely: 

• * 2G« 

Sij = 2G x eij - n Sj_j 
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ZGiKe 3 sinh 




2G z (eij 




CEq. 41 


This expression is intended to be used to predict the behavior of thermoplastic 
matrices in composites. The physical interpretation of the various parameters can be 
recapitulated as follows for the Eyring stress activated dashpot: 

K is the parameter which accounts for the activation energy; energy needed 
to overcome potential energy barrier. 

a is the parameter which accounts for the activation volume; the volume of 
the polymer segment which has to move as a whole in order for flow to 
take place . 

0 is the parameter which describes the pressure effect on shear 'yield* 

stress and determines the difference between tension and compression. 
B is the elastic bulk modulus 
Gj. is the shear modulus before 'yield* 

G 2 is the shear modulus after 'yield* 
n is the solid state viscosity 

Combining the hydrostatic and deviatoric components of the model, after a 
certain amount of algebra, an expression for the total strain rate in terms of 
total stresses was obtained. 

This is 


i = 
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Implementation in a Finite Element Code. 

The ABAQUS finite element code C211 which includes the option of a ''User 
Defined Subroutine - UMAT" , was selected as the platform for implementing the 
constitutive equation . The basic problem is to find the state of equilibrium 
corresponding to applied loads in which the resultant forces must vanish: ie, 
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{R(t)}-{F(t)J = O 

Where (R(t) } represents the forces due to externally applied loads 
{F(t) } represents the forces due to internal stresses 
The actual analysis process requires iteration, and in the general case, increments 
taken in terms of load or time steps. 

For a given increment, the main program assembles the global system of 
equations , solves for the unknown increment in displacements , evaluates the 
corresponding increment in strains and calls UMAT for each material point examined, 
to first calculate the increment in stresses according to the user's law, and then 
evaluate the material stiffness matrix. At this point, the program passes into UMAT 
the variables that need to be defined and/or updated and those needed for the 
calculations . 

The user must update the stress vector at the material point in accordance with 
the constitutive law. Another set of variables that must be defined for use in the 
next iteration is the material Jacobian relation between the increment in stresses and 
strains. This is used by ABAQUS to evaluate the element stiffness matrices in the 
mesh. 

The coding of UMAT is based on the following logic: 

In a given time step, dt=t(n) -t(n-l) : 

(1) From ABAQUS we know: 

the stresses {cr) and (S) 
the strains {6} and {e} at t=t(n-l) 
and an estimate of the increment in strains {a«} corresponding to dt. 

(2) The constitutive equation is expressed in terms of {S} and { e] . 

(3) The trapezoidal finite difference approximation is used for time integration 
to obtain a new relation between (AS) and {as}. 

(4) The six nonlinear equations for {AS} are solved using Newton -Raphson iteration. 
(The resulting algebraic equations are solved separately and not as a system, for 
reasons to be mentioned later ) . 

(5) The stress vectors {S } and {<r} are updated. 

3 act 

(6) The Jacobian matrix C~ — 3 is then defined. 

d AC 

The details of the manipulations required for this process are extensive C223 and are 
not presented here. 

Results Using Examples . 

The implementation of the time -dependent constitutive law requires careful 
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handling. To begin a set, preliminary tests were made for a simple 1-D linear 
viscoelastic model, to check the behavior of the subroutine as the applied loading rate 
changes and parameters become nonlinear. The manner in which variables pass into 
and out of UMAT under these circumstances was of particular interest. Specific 
results are not presented here; the experience gained during these tests, however, 
they did indicate the possibility of instabilities in the time integration procedure if the 
time steps are not carefully monitored. 

Predictions using a uniaxial version of the model and a simple numerical analysis 
program were then obtained. Two different sets of material properties were used to 
illustrate the sensitivity of the model to some of them. (See Table II-E-l) 

Table n-E-1 
MATERIAL PROPERTIES 





SET 1 

SET 2 

1 

K 

i 

5.0E-18 | 

1.0E-20 

t 

a 

i 

9.0E-08 | 

9.0E-08' 

1 

0 

i 

1 * OE - 14 | 

2.0E-08 

1 

Gl 

i 

. 88E+10 | 

* 88E+10 

1 

G2 

i 

. 88E+09 1 

• 88E+09 

1 

B 

i 

4.1E+10 1 

4 ♦ 1E+ 10 

1 

n 

i 

l.OE+16 ! 

1 . OE+16 

UNITS: 

dynes 

- cm - sec 




ldyne/cm* = O.lPa 



Theory and experiments were devised to show the influence of 

(1) Proper rate -dependence on 'yielding' 

(2) 'Yield' stress in compression > 'yield' stress in tension (by 'yielding', 
here , we mean the onset of nonlinearity . ) 

The first set of results shows that for three different strain rates, one obtains 
three different response -curves with common elastic branch, different 'yield' -stress 
level and almost parallel post -’yield' behavior; i.e. for higher strain -rate we get 
higher response. The comparison of these curves with experimental data is reasonably 
good (Fig. n-E-1). Similar results are obtained for uniaxial tension vs. 
compression (Fig. n-E-2). The behavior is also altered by adding constraints (Fig. 
n-E-3) . The last two cases illustrate the effect of the hydrostatic component on the 
overall behavior. 

The remaining examples are from the ABAQUS implementation of the constitutive 
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equation . The first example problem is a sample of a homogeneous , isotropic 
thermoplastic matrix ( Polycarbonate -Lexan ) subject to plane -strain uniaxial tension and 
compression via applied prescribed displacements at the one edge and proper boundary 
conditions at the other (shown below). The fibers are not included yet in the 
analysis, for reasons mentioned previously. A finite element mesh of 8-noded 
quadrilaterals with a 3x3/2x2 integration scheme was chosen from the ABAQUS element 
library to model the sample. 

Examining the behavior of the model in cycling loading (Fig. II -E -4), one can 
see the so-called 'ratcheting effect', i.e. the progressive compression of the 
hysteresis loops resulting from unloading in the post- 'yield' area. The unloading 
curve becomes almost immediately elastic, ie parallel to the original ascending branch. 
The new loading curve then follows the same path , until it reaches the onset of 
nonlinearity which is located at almost the original 'yield' stress level. 

When performing a couple of full cycles (tension and compression) under strain - 
control (Fig. II -E -5), no 'ratcheting' was allowed and despite the different stress 
level in opposite sign loadings, no "Baushinger effect" was noticed. After the first 
cycle was totally reversed, it followed exactly the same path 'as it did originally. 

Several load -and -hold tests were performed for uniaxial tension. We loaded very 
quickly up to a specific stress level which was on the elastic branch of the response 
or a little bit above the onset of ’yielding’, and then kept this stress constant through 
a certain time period. The plots of strain vs. time (Fig. II -E -6), show the 
expected creep behavior when the stress level is beyond the onset of nonlinearity. 

The final example problem was a beam bending test under pure moment . The 
stress distribution at the center cross-section of the beam was monitored as the value 
of the applied moment increased . When the tensile stress reached the onset of 
nonlinearity, the original linear stress distribution became curvilinear in the part of 
the beam which was under tension. This caused the over -shifting of the neutral axis 
towards the compression part, because now, more of the beam depth must be under 
tension, to satisfy the cross-sectional force and moment equilibrium equations: 

f h 

2 o^dy = 0 

ill 
z 

rh 

2 CT x y dy = 0 

-h 


and 


stress rate=1 MPa/sec 
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Figure II-E-5 

ress-Strain Relations in Two Full Cycles Under Strain Control. 
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When the stress overpasses the compressive ’yielding' point, its distribution over 
the beam depth becomes non-linear in both the outer tensile and compressive areas. 
The progressive overshifting of the neutral axis plotted versus the applied moment and 
some characteristic stress diagrams are shown in Fig. II-E-7. 

The Periodic Hexagonal Array (PHA) Mixing Model for Thermoplastics 

Several sophisticated mixing models have been developed recently in order to 
model, in a realistic way, the inelastic behavior of fibrous composites. The basic 
argument for using mixing models in composites is the level of scale that the user 
wants to look at. For small scale systems, where fiber -matrix interaction, local 
plastic deformations, debonding, delamination or even random fiber distribution cannot 
be ignored, use of the smoothed properties of a mixing model is inappropriate. On 
the other hand , when the macroscopic behavior of a composite system is under 
consideration, a good mixing model is appropriate. The combination of both scale 
levels in the same analysis is a difficult task, which requires special attention in the 
effort to capture those phenomena which affect the overall behavior. We belive, 
however, that the use of our constitutive relation for the ptatrix, within a reliable 
mixing model , will provide a useful tool in those areas where only macroscopic 
response is needed. 

The Periodic Hexagonal Array (PHA) model developed by G. J. Dvorak and J. 
L. Teply C233 for the elastoplastic behavior of metal matrix composites was selected 
for this work. The nonlinear matrix behavior is to be predicted by the constitutive 
relation presented in the preceeding subsections. 

The structure of the PHA is suitable for Finite Element Analysis, because it can 
be very easily used as a User Defined Subroutine to provide the characteristics of the 
material behavior in a general purpose F.E. program. Figure II -E -8 illustrates the 
interface between the F.E. program (ABAQUS) , the subroutine which defines the 
material behavior (UMAT), the mixing model (PHA) and the matrix constitutive relation 
(VISCOUS) . 

Our efforts are focused on the evaluation of the instantaneous material Jacobian 
matrix using the VISCOUS matrix constitutive relationship. Once this is determined we 
can estimate the increment in stresses by a simple multiplication with the given 
increment in strains. This material matrix, corresponding to smoothed properties, is 
obtained from an analysis of the representative volume element ( RVE ) , which requires 
a local finite element analysis of the mesh (fiber and matrix subelements) in the 
(RVE), subjected to boundary conditions which exclude rigid body motion [241. 



Shift in Neutral Axis (cm) 



Compressive Stress 
Tensile Stress 


F.E. program 

ABAQUS 
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The detailed steps required are as follows: 

(1) Given the overall stress, strain and strain increment, calculate local stress, 
strain and strain increment for each of the fiber and matrix subelements in the RVE. 
This is accomplished by use of stress and strain concentration factors . Their 
evaluation is presented, in detail, in [241. In fact, this step is performed only for 
matrix subelements, since the fiber behavior is assumed to be linear elastic and its 
material matrix does not depend on the stress -strain state. 

(2) Evaluate the fiber contribution to stiffness, using the elastic material matrix and 
the shape functions for the fiber subelements. This step is performed only once, 
since the above contribution remains constant. 

(3) Evaluate the matrix contribution to stiffness, using subroutine VISCOUS* for the 
material matrix and the shape functions for the matrix subelements in the RVE. This 
step has to be performed at every time step and globeil iteration, since the matrix 
behavior depends on the current stress -strain state . 

(4) Add the contribution to both fiber and matrix and come up with an expression of 
the overall instantaneous stiffness in terms of local moduli and volume fraction. 

(5) Multiply the overall instantaneous stiffness by the given increment in strains to get 
the corresponding increment in stresses. Return those values back to ABAQUS main. 

The fact that the proposed constitutive relation is defined in strain -space (given 
strain - return stress) makes it easy to incorporate into the PHA, since it is totally 
compatible with the displacement -based finite element method. 


* The subroutine VISCOUS is called at each integration point of the structure, for 
each one of the matrix subelements in RVE. It returns the material matrix 
(stress -strain relation), given the local stress -strain state and the current strain 
increment. The algorithms included in VISCOUS are exactly those described earlier in 
this section, representing the nonlinear time -dependent constitutive relation for 
thermoplastics . 
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PLANS FOR UPCOMING PERIOD 

Efforts during the upcoming period will be devoted to continuing the development of 
numerical analysis capabilities and applying them to specific experimental 
configurations. This work will be focused on the following areas : 

1) Completing the PH A model for thermoplastic composites. 

2) Investigating the effect of each parameter and stress/strain component and its 
significance on overall composite material behavior. 

3) Using a more sophisticated time integration operator in the model formulation, to 
insure efficiency, accuracy and rapid convergence. 

4) Developing a 2-D version of the PHA mixing model, in order to avoid performing 
3-D analyses. This task may involve considerable work if we reformulate the model, 
rather than just enforcing the constraints of the plane cases to the more general 3-D 
formulation . 

We have just begun investigating the problem of fracture toughness of 
thermoplastic composites by modeling Mode I and Mode II fracture specimens and 
determining the local stress fields around the crack tip. The interaction between 
fibers and thin matrix films close to the crack tip is expected to be significant, in 
particular when using a realistic time -dependent model for the matrix. Of course, the 
effort to combine a macroscopic analysis (say the Double cantilever Beam (DCB) test) 
and a close look at the crack tip , will not be easy . Comparisons between the 
classical fracture mechanics approach based on am infinitely sharp crack, and the kind 
of blunt cracks which appear to be encountered with thermoplastics, will also be made. 
Energy considerations will be included to determine the amount of energy that is 
dissipated versus that which is stored. Consideration will also be given to the bulk 
(hydrostatic) energy component versus the shear (deviatoric) . The final goal is to 
understand why thermoplastic composites are structurally tougher than those with 
conventional thermoset matrix materials. 

These analyses are being closely coordinated with experiments being carried out 
by Professor Stemstein. 

PUBLICATIONS AND PRESENTATIONS BY PROFESSOR M. SHEPHARD ON THIS SUBJECT 

"The N -Criterion for Predicting Crack Growth Increment", with N.A.B. Yehia, Engng. 
Fracture Mech. , Vol. 26, No. 3, 1987. 
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•‘Progress on Automated Finite Element Modeling", presentation at Pratt and Whitney 
Aircraft, Hartford , CT, September 18, 1986. 

"Automated Finite Element Modeling", Grumman Aircraft, Bethpage, MY, February 6, 
1987. 

"Nonlinear Finite Element Modeling of Composites”, ONR review of SDI related 
composites work, U. of Maryland , College Park, Maryland , March 31, 1987. 

"Finite Element Modeling of Thermoplastic Composites - Constitutive Relations for 
Matrix Materials”, with R.J. Bankert, S.S. stemstein and N.D. Lambropoulos , 
abstract accepted for ASTM Symposium on Advances in Thermo -plastic Matrix Composite 
Materials”, October 1987. 
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P. GENERAL BEAM THEORY FOR COMPOSITE STRUCTURES 


Sr. Investigator: O. Bauchau 


INTRODUCTION 

This work has concentrated on the development of a general beam theory for thin 
walled structures , including cross -sectional deformations. Both numerical and 

experimental aspects have been investigated. Since no assumption is made in this 
work about the cross-sectional deformation, the resulting analysis should yield results 
equivalent to a full three-dimensional model of the structure. The cost of the 
analysis, on the other hand, should be orders of magnitude lower them a conventional 
three dimensional finite element analysis. 

STATUS 

The proposed generalized beam theory has important applications in light weight 
aeronautical structures . Wing boxes and fuselages are often composed of flat or 
curved panels, reinforced by members with fairly deep stiffness. A major problem in 
the design of such built-up Shell structures subjected to in -plane compressive and/or 
shearing loads is the existence of many different failure modes, and the possible 
interaction between these modes. The buckling and poet -buckling analysis of these 
structures using the finite element method has become popular, but the non-linear 
analysis of complex shell structures requires very large amounts of computing time and 
is too expensive to be used systematically in the predesign process. The increasingly 
widespread use of composites for thin -walled structures adds impetus to the effort to 
develop efficient methodologies which make no a priori assumptions regarding failure 
modes; particularly because of the relatively limited experience with such applications. 

PROGRESS DURING REPORTING PERIOD 

During this research period, the generalized beam model was developed. The 
governing equations are as follows: 

-Ku" + Ku* + Ku = Q Eq. (1) 

Where: 

K, K, and K are stiffness matrices resulting from the 

discretization of the cross-section of the structure, 

u(z) is the vector of modal displacements, 

( )* denotes a derivative with respect to the axial variable, z, 
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Q is the load vector. 

Equation (1) represents a set of coupled differential equations that can be solved 
by first finding the eigen deformation modes from the following quadratic eigenproblem . 

( -p£k + PjK + K) Ui = 0 Eq. (2) 

Where: 

Ui are the eigen deformation modes which characterize the in-plane 
deformation modes of the section, and 

Pi are the associated eigenvalues. 

The task of solving Eq. (2) proved to be extremely difficult: a special 

procedure based on the Generalized Lanczos Algorithm was developed, and the eigen 
deformation modes for various sections were successfully obtained. The convergence 
of the eigenvalue extraction routine , however , was fairly slow , not because of an 
inherent deficiency of the algorithm, but because of the existance of closely -spaced 
eigenvalues . This means that the solution of the overall problem , based on the 
superposition of these modes , would be very slow in converging as well , since the 
contribution of each mode is inversely proportional to the valdes of the eigennumbers . 

Another approach to the solution of Eq . ( 1 ) has now been taken . The 

displacement function vector u is expanded in a series of Tchebychev polynomials 

along the axis of the beam. This numerical procedure is showing great efficiency. 

PLANS FOR THE UPCOMING PERIOD 

Following guidance from NASA/AFOSR monitors, on the occasion of the program 
site visit on December 18 and 19, 1986, support for this project under the subject 

grant was terminated as of the end of the reporting period. In view of the promise 

of the Tchebychev polynominal approach, however, the methodology development for 
the efficient analysis of composite structures with thin cross -sections is expected to 
be completed under other funding. 


PART m 

TECHNICAL INTERCHANGE 
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TECHNICAL INTERCHANGE 


Technical meetings, both on- and off-campus, provide for the interchange of 
technical information. In order to assure that Rensselaer faculty /staff members can 
participate , a central listing of upcoming meetings is compiled , maintained and 
distributed on a periodioc basis. The calendar for this reporting period is shown in 
Table III-l. Table HI -2 shows the meetings attended by RPI composites program 
faculty /staff /students during the reporting period. Some on -campus meetings, with 
special speakers particularly relevant to composites, are listed in Table HI -3. A list 
of composites -related visits to relevant organizations, attended by RPI 
faculty/staff /students , along with the purpose of each visit is presented in 
Table m-4. 

The diversity of the research conducted within this program has increased over 
the last several years. To insure information transfer, once -a -week luncheon 
programs have been held among the faculty and graduate students involved (listed in 
Part IV. Personnel - of this report). These meetings are held during the academic 
year and are known as "Brown Bag Lunches" (BBL's) , since attendees bring their own. 
Each BBL allows an opportunity for graduate students and faculty to briefly present 
plans for, problems encountered in and recent results from their individual projects. 
These seminars also are occasions for short reports on the content of off-campus 
meetings attended by any of the faculty /staff participants (see Tables III -2 and 
III-4) and for brief administrative reports, usually on the part of one of the 
Co -Principal Investigators. Off-campus visitors, at RPI during a BBL day, are often 
invited to "sit in". Table IH -5 lists a calendar of internal, oral progress reports as 
they were given at BBL's during this reporting period. 

As an important part of the steps taken to increase communication between NASA 
researchers and their RPI counterparts in areas of interest under this grant, a series 
of Research Coordination Meetings have been held with members of NASA Langley and 
Lewis Research Center's Materials and Structures scientist/engineers. These meetings 
are summarized in Table III -6. 


* 
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I 

Table III-l 

COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
Calendar of Composites -Related Events 


May 1, 1986 through April 30, 1987 


DATES 

MEETING 

SPONSOR 

PLACE 

12-19 May 
86 

workshop on Failure Mechanics 

ONR 

College Park, MD 

14 May 86 

Symposium on Advanced 
Composite Materials 

GE 

Schenectady, NY 

19-21 May 
86 

27th Structures, Structural 
Dynamics & Materials Conf. 

AIAA/ASME/ 

ASCE/AHS 

San Antonio, TX 

19-23 May 
86 

Inti. Conf. on Advanced 
Composite Materials 

Chinese Soc. 
for Matl. Sc. 

Taipei, Taiwan 

21 May 86 

An Introduction to 
Advanced Compos ites - 
Materials Processes, 
Equipment & Applications 

CoGSME 

Philadelphia, PA 

22 May 86 

Advanced RTM for Advanced 
Composite Production 

CoGSME 

Philadelphia, PA 

2-3 Jun 
86 

Compos ites Consortium 
Program Review 

SDIO/ONR 

Woods Hole, MA 

2-6 Jun 
86 

Inti. Conf. on Role of 
Fracture Mechanics in 
Modern Technology 

Kyushu Univ. 

Fukuoka, Japan 

4-5 Jun 
86 

Technology in 1990* s: 
Advanced Materials 

SME 

London, England 

8-13 Jun 
86 

Spring Conference on 
Experimental Mechanics 

SEM(SESA) 

New Orleans , LA 

9-12 Jun 
86 

AUT0C0M *86 - Advanced 
Applications of Composites 
for Automotive 

SME 

Dearborn , MI 

11-13 Jun 
86 

3rd Symposium on Nonlinear 
Constitutive Relations for 
High Temperature Applications 


Akron , OH 

16-20 Jun 
86 

10th U.S. National 
Congress on Applied 
Mechanics 

ASME 

Austin, TX 
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Table III-l (continued) 

COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
Calendar of Composites -Related Events 

May 1, 1986 through April 30, 1987 


DATES MEETING SPONSOR PLACE 


23-25 

Jun 86 

3rd Japan-U.S. Conference 
on Composite Materials 

NSF 

Tokyo, Japan 

30 Jun - 
2 Jul 86 

19th National Symposium 
on Fracture Mechanics 

ASTM 

San Antonio, TX 

30 Jun - 
4 Jul 86 

4th Inter. Carbon Conf. 

D.K.G. 

Baden-Baden, 
west Germany 

21-23 Jul 
86 

2nd Inti. Symposium on 
the Nondestructive 
Characterization of 
Materials 

Industrial 
Matls. Res. 
Inst. (NRCC) 

Montreal, Canada 

19-21 Aug 
86 

Conf. on Nondestructive 
Testing & Evaluation of 
Advanced Matls a Composites 

DoD 

Colorado Springs, 
CO 

25-29 Aug 
86 

Inti. Conf. & Exposition 
on Engineering Ceramics 

Am. Soc. for 
Metals 

Buffalo, NY 

25-29 Aug 
86 

COMP '86 Symposium - 
Engineering Applications 
of New Composites 

Univ. of 
Patras 

Patras, Greece 

1-5 Sep 
86 

Symposium 

IUTAM 

Paris, France 

7-12 Sep 
86 

ACS 192nd Natl. Mtg. 

ACS 

Anaheim, CA 

8-11 Sep 
86 

Fabricating Composites 
*86 Conf. & Exposition 

CoGSME & SME 

Baltimore, MD 

22-26 Sep 
86 

1st World Congress on 
Computational Mechanics 
(WCCM) of the Inti. 
Assoc, for Computational 
Mechanics (IACM) 

Univ. of TX 

Austin, Tx 

24 Sep 
86 

15th Videoconference 
In-Process Control for 
Manufacturing 

IEEE Student 
Branch, RPI 
ECSE/CMPTT 

Troy, NY 
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Table III-l (continued) 

COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
Calendar of Composites -Related Events 

May 1, 1986 through April 30, 1987 


DATES MEETING SPONSOR PLACE 


4-5 Oct 
86 

1986 Materials Science 
Seminar "Computer 
Simulations in Materials 
Science" 

ASM MSD 

Orlando, FL 

4-9 Oct 
86 

ASM Materials Week 

ASM 

Lake Buena Vista, 
FL 

7-9 OCt 
86 

1st Conf. on Composite 
Materials 

ACS ft Univ. 
of Dayton 

Dayton, OH 

7-9 OCt 
86 

18th Natl. SAMPE Tech. 
Conf.: Matls for Space: 
The Gathering Momentum 

SAMPE 

Seattle, WA 

8-10 Oct 
86 

0RCAL '86 (Orange Cty 
Manufacturing a Metal- 
working Conf. & Expo.) 

ASM & SME 

Anaheim, CA 

10 OCt 
86 

Fastening Advanced 
Composites Conf. 

SME 

Renton, WA 

14-16 OCt 
86 

18th National 
Technical Conf. 

SAMPE 

Seattle, WA 

20-22 OCt 
86 

Aircraft Systems, 
Design & Tech. Mtg. 

AIAA/AHS/ASEE 

Dayton , OH 

28-30 OCt 
86 

MMC Spacecraft 
Survivability 
Tech. Conf. 

DoD 

Menlo Park, CA 

2-5 NOV 
86 

Optical Methods & 
Composites 

SEM 

Keystone, 00 

3 Nov 86 

Test Methods ft Design 
Allowables for Fiber 
Composites: 2nd Symposium 

ASTM 

Phoenix, AZ 

10 Nov 86 

Southwest Mechanics 
Lecture Series 

♦ 

Norman, OK 

10-14 Nov 
86 

Inti. Conf. ft Expo, on 
Castings 

ASM 

Chicago, IL 
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Table Ill-l ( continued ) 


COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
Calendar of Composites -Related Events 

May 1, 1986 through April 30, 1987 


DATES 

MEETING 

SPONSOR 

PLACE 

2-4 Dec 
86 

Composites Materials: 
Analysis Testing ft Design 

SEM 

Chicago, IL 

7-12 Dec 
86 

Winter Annual Mtg. 

ASME 

Anaheim, CA 

7 Jan 87 

Inti. Conf. on Constitutive 
Laws for Engineering Materials: 
Theory and Applications 

Tucson, AZ 

12-16 Jan 
87 

Gordon Conference 
on Composites 

Gordon 

Conf. 

Santa Barbara, CA 

19-22 Jan 
87 

Composites in Manufacturing, 
6th Conference ft Exposition 

CoGSME 

Anaheim, CA 

11 Mar 87 

workshop on Composite 
Materials - Interface 
Science 

ONR 

Leesburg, VA 

30-31 Mar 
87 

ONR Review of SDI 
related composites 
work 

ONR 

College Park, MD 

6-8 Apr 
87 

28th Structures, Struc- 
tural Dynamics ft 
Materials Conf. 

AIAA/ASME/ 

ASCE/AHS 

Monterey, CA 

6-9 Apr 
87 

32nd Inti. SMPE Symposium 
ft Exhibition 

Society of 
Advanced 
Matls ft 
Processing 
Engineering 

Anaheim, CA 

8 Apr 
87 

Composites : The Future 

is Now 

Capital 
Region Tech. 
Development 
Council 

Troy, NY 

9-10 Apr 
87 

Dynamics Specifications 
Conference 

AIAA 

Monterey, CA 

10 Apr 
87 

Supportability Seminar 

CoGSME 

Anaheim, CA 
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DATES 
12-13 
May 86 


14 May 86 


19-21 May 
86 


2-3 Jun 
86 


2-4 Jun 
86 

11-13 
Jun 86 


Table III -2 

COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
Pertinent Professional Meetings Attended 

May 1, 1986 through April 30, 1987 


MEETING 

(MIR Workshop on Failure Mechanics (Prof. Dvorak), College 
Park, MD 

Prof. Dvorak presented the paper: 

“Fracture Mechanics of Metal Matrix Composites” 

Symposium on Advanced Composite Materials (Prof. Dvorak), 
Schenectady, NY 

Prof. Dvorak presented the paper: 

“Analysis of Fatigue Cracking of Fibrous Metal Matrix 
Laminates” 

27th Structures, Dynamics & Materials Conference (Prof. Loewy), 
San Antonio, TX 

Prof. Loewy presented the paper: 

“Application of the Principal Curvature Transformation to 
Nonlinear Rotor Blade Analysis" 

SDIO/ONR Composites Consortium Program Review (Prof. Dvorak), 
Woods Hole, MA 

Prof. Dvorak presented the paper: 

“Damage in Metal Matrix Composites” 

AHS Annual Forum (Prof. Loewy), Washington, DC 


3rd Symposium on Nonlinear Constitutive Relations for High 
Temperature Applications (Prof. Krempl), Akron, OH 

Professor Krempl presented the papers: 

“The Viscoplasticity Theory Based on Overstress Applied to 
the Modeling of Nickel Base Super Alloy at 815°C“ 

“Cyclic Uniaxial and Biaxial Hardening of Type 304 Stainless 
Steel Modeled by the Viscoplasticity Theory Based on 
Overstress” 

“A Simplified Orthotropic Formulation of the Viscoplasticity 
Theory Based on Overstress" 
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DATES 
16-20 
Jun 86 


23-25 
Jun 86 


30 Jun - 
4 Jul 86 


1-5 sep 

86 


21-22 Oct 
86 


10 Nov 86 


7-12 
Dec 86 


Table in -2 (continued) 

COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
Pertinent Professional Meetings Attended 

May 1, 1986 through April 30, 1987 


MEETING 

10th U.S. National Congress of Applied Mechanics (Prof. 
Dvorak), Austin, TX 

Prof. Dvorak presented the paper: 

"Thermal Expansion of Elastic -Plastic Composite Materials" 

3rd Japan-U.S. Conference on Composite Materials (Prof. 
Diefendorf), Tokyo, Japan. 

Professor Diefendorf presented the paper: 

"The Relationship of Structure to Properties in Carbon Fibers" 

4th International Carbon Conference (Prof. Diefendorf), 
Baden-Baden, Deutschen Keramischen Gesellschaft, West Germany 

Professor Diefendorf presented the paper: 

"The Chemical Vapor Deposition of Carbon Capillary Tubes II” 

IUTAM Symposium (Prof. Dvorak), Paris, France 

Prof. Dvorak gave the Invited Lectures: 

"Thermomechanical Couplings in Solids” 

"Thermomechanical Deformation and Coupling in Elastic-Plastic 
Composite Materials" 

aiaa/ahs/asee Aircraft Design Mtg (Prof. Loewy), Dayton, OH 
Prof. Loewy was a Panel Discussion Member 

Southwest Mechanics Lecture Series (Prof. Krempl), Norman, OK 
Professor Krempl presented the paper: 

"Biaxial Fatigue and Deformation Behavior of Graphite/Epoxy 
Composites" 

ASME Winter Annual Meeting (Prof. Dvorak), Anaheim, CA 

Prof. Dvorak gave the Invited Lecture: 

"Fatigue Damage Analysis in Metal Matrix Laminates" 
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Table III -2 (continued) 


COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
Pertinent Professional Meetings Attended 

May 1, 1986 through April 30, 1987 


DATES 
7 Jan 87 


12-16 Jan 
87 


11 Mar 87 


30 Mar 87 


31 Mar 87 


6-7 Apr 
87 


MEETING 

International Conference on Constitutive Laws for Engineering 
Materials: Theory and Applications (Prof. Krempl), Tucson, AZ 

Prof. Krempl presented the paper: 

"Isotropic and Orthotropic Formulations of the Viscoplastic 
Theory Based on Overstress" 

Gordon Conference on Composites (Prof. Sternstein), Santa 
Barbara, CA 

Professor Sternstein gave the Lecture: 

"Thermoplastic Matrix Composites" 

ONR Workshop on Composite Materials - Interface Science 
(Prof. Dvorak), Leesburg, VA 

Prof. Dvorak presented the paper: 

"Cracks Approaching Interfaces: The Image Crack Method” 

ONR SDI Review (Prof. Dvorak), College Park, MD 
Prof. Dvorak presented the paper: 

"Analysis of Metal Matrix Composites for Spacecraft 
Applications” 

ONR Review of SDI Related Composites Work (Prof. M. Shephard), 
College Park, MD 

Professor Shephard gave the Lecture: 

"Nonlinear Finite Element Modeling of Composites” 

28th Structures, Dynamics & Materials Conference 
(Prof. Loewy), Monterey, CA 



65 


Table m-3 

COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
Composites -Related Meetings/Talks Held at RPI 

May 1, 1986 through April 30 > 1987 


SUBJECT SPEAKERS (RPI) DATE 

SHORT COURSE: Prof. O. Bauchau 7/21-25/86 

Advanced Composite Prof. R. Diefendorf 

Materials and 
Structures 

Structure Failures ft L. Coffin 9/2/86 

Their Impact on 
Engineering Knowledge 


Workshop on Composite 
Materials and 
Structures for 
Rotorcraft 


inhomogeneous Swelling 
ft Sol vents/ So lids/S tress 
interactions in Multiphase 

S . Sternstein 

10/2/86 

Damage Monitoring, Life 
Prediction ft Life Ex- 

L. Coffin 

10/21/86 

tension of Engineering 
Structures 

SPEAKERS 


NASA/AFOSR site Visit 

0. Bauchau 

12 / 18 - 19/86 

and Program Review 

R. Diefendorf 
G. Dvorak 
E. Krempl 

R. Loewy 
V. Paedelt 
M. Shephard 

S. Sternstein 
S. Winckler 



o. Bauchau 
M. Da r low 
R. Diefendorf 
R. Loewy 
s. Winckler 


9/10-11/86 
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Table HI- 3 (continued) 

COMPOSITE MATERIALS AMD STRUCTURES PROGRAM 
Composites -Related Meetings/Talks Held at RPI 

May 1, 1986 through April 30, 1987 


SUBJECT 
Theory of 
Vi scopl as t ic i ty 
Based on Over- 
stress with 
Applications 

Minimizing 
Residual Stresses 
in Injection 
Molded Parts 

Pultruded 
Fiberglass -A 
Structural 
Alternative 

Failure Criterion 
of Fiber 
Reinforced 
Plastics and 
Optimum Fiber 
Orientations 

Servo -Control led 
Testing and 
Materials 
Modelling 

Localization 
of Plastic 
Deformation 

Presentation of 
U.S. Army Engine 
Developmental 
Research Program 


SPEAKERS 
David Yao 

Doctoral Dissertation 


Ming J. Liou 
M. I. T. 
Cambridge, MA 


Russ Whitfield 

Morrison Molder Fiber Glass Co. 
Bristol, VA 


Prof. K. Ikegami 

Tokyo institute of Technology 

Tokyo, Japan 


Dr . E . Krempl 
R. P. I. 


Prof. A. S. Douglas 
Johns Hopkins Univ. 


Dr. J. Acurio 


DATE 

2/11/87 


2/12/87 


2/27/87 


2/27/87 


3/23/87 


3/24/87 


4/15/87 
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Table III-4 


COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
Composites -Related Visits to Relevant Organizations 


Faculty Member 
G. Dvorak 


G. Dvorak 


G . Dvorak 


R. Diefendorf 


0. Bauchau 
R. Diefendorf 
R. Loewy 


R. Loewy 


M . Shephard 


E. Krempl 


G. Dvorak 
R. Diefendorf 

R. Loewy 

S. Sternstein 


May 1, 1986 through April 30, 1987 


Purpose of Visit 
Colloquium: "Damage 
Mechanics of Composite 
Materials” 

Presented paper: 
"Fracture Mechanics 
of Metal Matrix 
Composites” 

Symposium: "Analysis 
of Fatigue Cracking 
of Fibrous Metal 
Matrix Laminates” 

Presented Seminar: 
"Carbon Fibers From 
Mesophase Pitch” 

Tour of Facilities 
& Technical Discussion 


Discussion of 
Structural Dynamics 
Research 

Presented Seminar: 
"Progress on Automated 
Finite Element 
Modeling" 

Presented Seminar: 
"Biaxial Fatigue and 
Deformation Behavior 
of Graphite Epoxy 
Composites" 

Discussion of Langley 
R. C. research program 
in composites 


Location 

Northwestern Univ 
Evanston , IL 


U. of MD 

College Park, MD 


G.E. 

Schenectady , NY 


Rockwell Science 
Ctr. , Thousand Oaks, 
CA 

U.S. Composites 
Corp. , Rensselaer 
Technology Park, 

No . Greenbush , NY 

Jet Propulsion Lab. 
Pasadena, CA 


Pratt & Whitney 
Aircraft , 
Hartford, CT 


University of 
Delaware, 
Newark , DE 


NASA Langley R.C. 
Hampton, VA 


Date( s 1 
5/2/86 


5/12-13/ 

86 


5/14/86 


7/11/86 


7/29/86 


8/6/86 


9/18/86 


9/19/86 


11/3/86 



68 


Table III -4 (continued) 



COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
Composites -Related Visits to Relevant Organizations 


1987 

May 1, 1986 through April 30, 


1 

Faculty Member 
E. Krempl 

Purpose of Visit 
Presented Lecture; 
"Biaxial Fatigue and 
Deformation Behavior 
of Graphite/Epoxy 
Composites" 

Location 
University of 
Oklahoma, 
Norman, OK 

Datef s ) 
11/10/86 

G. Dvorak 

Colloquium: "Plasticity 
of Composite Materials" 

Texas A&M Univ 
College Sta, TX 

11/11/86 

G . Dvorak 

Colloquium: "Plasticity 
of Composite Materials" 

Rice University 
Houston, TX 

11/12/86 

G . Dvorak 

Colloquium: "Recent 
Developments in 
Plasticity of Fiber 
Metal Matrix Composites 

Yale Univ. 

New Haven , CT 

1/21/87 

G. Dvorak 

Symposium: "Bimodal 
Plasticity Theory of 
Composite Materials" 

U. of Florida 
Gainesville, FL 

1/28-30/ 

87 

M. Shephard 

Presented Seminar: 
"Automated Finite 
Element Modeling" 

Grumman Aircraft 
Bethpage , NY 

2/6/87 

G. Dvorak 

i 

1 

Presented Seminar: 
"Recent Developments 
in Plasticity of 
Composite Materials” 

U. of CA 
Berkeley, CA 

2/23/87 

S. Stemstein 

Presented Seminar: 
"Thermoplastic Matrix 
Composites" 

U. of Connecticut 
Storrs, CT 

3/8/87 

R. Diefendorf 

Presented Seminar: 
"Carbon Fibers” 

SUNY, Buffalo, NY 

4/2/87 

R. Diefendorf 

Presented Seminar: 
"Chemical Vapor 
Deposition" 

Sohio, 

Niagara Falls, NY 

4/6/87 

G. Dvorak 

Presented Seminar: 

"A Bimodal Plasticity 
Theory of Composite 
Materials” 

Brown Univ. 
Providence , RI 

4/6/87 


I 
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1987 

Faculty Member 
S. Sternstein 


R. Diefendorf 


Table III -4 (continued) 

COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
Composites -Related Visits to Relevant Organizations 

May 1, 1986 through April 30, 


Purpose of Visit 
Presented Seminars 
"Matrix Dominated 
Mechanical Properties 
of Composites" 

Presented Sach’s 
Memorial Lecture: 
"Composite Processing" 

Presented Seminar: 
"Recent Developments 
in Composites 
Plasticity" 


Location 

TRW Research Ctr 
Los Angeles, CA 


Syracuse University 
Syracuse, NY 


Lawrence Liver- 
more Labs . , 
Livermore , CA 


Date( s ) 
4/10/87 


4/21/87 


G. Dvorak. 


4/23/87 
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DATE 
02 -May 

05 -Sep 

12 -Sep 

19 -Sep 

26 -Sep 
03 -Oct 

10 -Oct 
17 -Oct 

24 -Oct 
31 -Oct 
07 -Nov 
14-Nov 


Table in -5 

COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
Brown Bag Lunch Schedule 

May 1, 1986 through April 20, 1987 


TOPIC 

Administrative Report 
Edge Failures 
Fabrication Technology 


RESP. FACULTY 

Loewy 

Sham 

Bundy/Hagerup/Paede 1 t 


Administrative Report 
Fracture Toughness in Thermoplastic 
Matrix Composites 
Failure in Metal Matrix Composites 


Loewy 

Sternstein 

Dvorak 


General Discussion - 
Rpt on RPI Workshop on Composites 
in Rotorcraft 


Loewy/Diefendorf/ 

Bauchau/Winckler 


Administrative Report 
Chemical Vapor Deposition 
Anisotropic Beam Theory 


Diefendorf 

Diefendorf 

Bauchau 


General Discussion - 

Preparation of Ceramic Fibers Interrante 


Administrative Report 

Time -Dependent Deformation of Comp. 

Numerical Analysis of Composites 


Loewy 

Krempl 

Shephard 


General Discussion - 

High Temperature Deformation Behavior Krempl 


Administrative Report 

Fabrication Technology 

Ceramic Matrices - Prep. & Properties 


Diefendorf 

Bundy/Hagerup/Paede 1 t 
Doremus 


General Discussion - 

Mechanics of Damage in Composites Sternstein 


Administrative Report Diefendorf 

Intermetallic Matrix Fabrication Tech. German 


General Discussion - 

Processing at Very High Temperatures Doremus 


Administrative Report 
Organometallic Precursors to Ceramics 


Loewy 

Interrante 
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DATE 
21 -Nov 

28 -Nov 

05 - Dec 

12 - Dec 
16 -Jan 

23 -Jan 

30 -Jan 

06 - Feb 

13 - Feb 

20 -Feb 
27 -Feb 

06 -Mar 
13 -Mar 


Table IH -5 (continued) 

COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
Brown Bag Lunch Schedule 

May 1, 1986 through April 20, 1987 


TOPIC 

General Discussion - 

Mechanics of Damage in Composites 

Thanksgiving Holiday (No Mtg) 

Administrative Report 
High Temperature Chemistry in 
Ceramic Matrix Composites 

General Discussion - 
Carbon/Carbon and CVD 

Administrative Report 
General Discussion: RP-3 

Design Evolution 

Administrative Report 

Anistropic Beam Theory 

Failure in Metal Matrix Composites 

General Discussion - 
Unsymmetric Laminates 

No Meeting 

Administrative Report 

Modal Analysis of RP-2 Sailplane 

Administrative Report 
Fabrication Technology 
Fracture Toughness in Thermoplastic 
Matrix Composites 

Administrative Report 
Guest Speaker 

Failure Criterion of Fiber Reinforced 
Plastics and Optimum Fiber 
Orientations 


RESP. FACULTY 
Krempl 


Diefendorf 

Wiedemeier 


Diefendorf 

Diefendorf 

Paedelt 


Loewy 

Bauchau 

Dvorak 


Winck 1 er/Bauchau 


Diefendorf 

Swaybill 

Loewy 

Paedelt 

Sternstein 


Loewy 

Prof. K. Ikegami 
Tokyo Institute 
of Technology 


Administrative Report 

Failure in Metal Matrix Composites 

Numerical Analysis of Composites 

General Discussion - 

Mechanics of Damage in Composites 


Diefendorf 

Dvorak 

Shephard/Lambropou los 


Sternstein/Krempl 


Table III -5 (continued) 


COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
Brown Bag Lunch Schedule 

May 1, 1986 through April 20, 1987 


DATE 
20 -Mar 

TOPIC 

No Meeting (Spring Break) 

RESP, FACULTY 

27 -Mar 

General Discussion - 

Project Progress on the U.R.I. Program 

Diefendorf 

03 -Apr 

Administrative Report 

Chemical Vapor Deposition 

Time -Dependent Deformation of Comp. 

Diefendorf 

Diefendorf 

Krempl 

10 -Apr 

Administrative Report 
Anistropic Beam Theory 
Progress in Failure Analysis 

Diefendorf 

Bauchau 

Sham 

17 -Apr 

General Discussion - 

High Temperature Deformation Behavior 

Krempl/Winckler 

24-Apr 

Administrative Report 
Fracture Toughness in Thermoplastic 
Matrix Composites 
Numerical Analysis of Composites 

Loewy 

Sternstein 

Shephard 


Visitor: Steven Smith, Courtaulds Ltd., England 

(Mfg’r. of Carbon Fibers) 
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Faculty Member 
R. Diefendorf 
G. Dvorak 

R. Loewy 

S. Sternstein 

R. Loewy 


R. Loewy 


R. Loewy 


R. Diefendorf 
G. Dvorak 

R. Loewy 

S . Sternstein 


Table III -6 

COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
Review of Research Center Interactions 

May 1, 1986 through April 30, 1987 


Purpose of Visit 
Discussion of Langley 
R. C. research program 
in composites 

Location 

NASA Langley R.C. 
Hampton, VA 

Date( s > 
11/3/86 

Discussion of Composites 
research with Dr. M. 
Greenfield 

NASA HDQS 
Washington , DC 

11/19/86 

Discussion of Composites 
research with Drs. A. 
Amos, G. Haritos and 
M. Salkind 

AFOSR HDQS 
Washington, DC 

11/19/86 

Discussion of Conqposites 
research with Dr. D. 
Mulville 

NASA HDQS 
Washington, DC 

1/12/87 

Discussion of Composites 
research programs at 
Lewis R.C. 

NASA Lewis R.C. 
Cleveland, OH 

2/6/87 
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During the week of July 21-25, 1966, RPI offered, for the seventh time, a 

special short course in composite materials and structures. Seventeen graduate 
engineers from government and industry enrolled. In addition to the RPI speakers 
shown in Table III-3, Prof. L. Phoenix of Cornell University, Mr. Bob Riley of 
McDonnell Aircraft and Dr . Stephen Tsai of the Air Force Material Laboratories 
lectured . The announcement brochure , listing lecturers and the subject matter is 

attached as an appendix to this report, and the participants and their organizations 
are listed in Table III- 7. 

The international workshop on composite materials and structures for rotor craft, 
also shown in Table III-3, was conducted at RPI on 9/10-11/86 at the suggestion of 
the Army Research Office. A list of attendees, including 61 representatives of 

industry, government and academia, is shown in Table III -8 and the agenda for the 
workshop is appended. This by -invitation -only meeting was considered sufficiently 
successful that follow-on workshops are being considered for future years. 


Table III -7 


COMPOSITE MATERIALS AND STRUCTURES PROGRAM 

Short Course: Composite Materials and structures 

Participants and Affiliations 

July 21, 1986 through July 25, 1986 


Wayne Burgess 
Design Engineer 
General Electric Co. 

Lakeside Avenue Rm. 1311 
Burlington, VT 05402 

Robert Filler 
Member Technical Staff 
Structures Dept. 

McDonnell Douglas Helicopter Co. 

4645 South Ash Avenue - H/S G3 
Tempe , AZ 85282 

Terry Gossard 
Associate Engineer 
Atlantic Research Corp. 

7511 Wellington Road 
Gainesville, VA 22065 

James K. Gran 
Research Engineer 
SRI International 
333 Ravenwoods Avenue 
Menlo Park, CA 94025 

Carl Husen 

Sr. Engineering Specialist 

Ford Aerospace & Communication Corp. 

Ford Road 

Newport Beach, CA 92658-9983 
Donald Johnson 

Engineering Structures Technician 
Boeing Vertol Co. 

P 24-03, P.O. Box 16858 
Philadelphia, PA 19142 

Dan Malwitz 
Engineer 

Contraves Goerz Corp. 

610 Epsilon Drive 
Pittsburgh, PA 15238 


Brian McKillop 
Engineer 

Fiber Materials, Inc. 

Biddeford Industrial Park 
Biddeford , ME 04005 

Steve Nelson 

Stress Engineer 

Heath Tecna Aerospace Co. 

19819 84th Avenue S. 

Kent , WA 98032 

Eric Preissner 
Aerospace Engineer 
Naval Air Development Center 
Code 6055 

Warminster, PA 18974 

Eric Rahnenfuehrer 

Material & Processing Engineer 

Kaman Aerospace Corp. 

P. 0. Box 2 

Bloomfield, CT 06002-0002 

Rebecca Sidler 
Engineer 

Chamberlain Mfg. Corp. 

550 Esther Street 
Waterloo, IA 60704 

Richard Soloski 
Member Technical Staff 
Structures Dept. 

McDonnell Douglas Helicopter Co. 
4645 South Ash Avenue - M/S G3 
Tempe, AZ 85282 

Allan Stanley 

Flight Test Development Engineer 
The 6420 Test Group - ENMT 
Edwards AFB, CA 93523 


Table m-7 (continued) 


COMPOSITE MATERIALS AND STRUCTURES PROGRAM 

Short Course: Composite Materials and Structures 

Participants and Affiliations 

July 21, 1986 through July 25, 1986 


Jordon Watters 
Engineer 

Lawrence Livermore National Labs. 
P. 0. Box 808 
Livermore, CA 94550 

Paul Wren 
District Manager 
TRW 

1400 Talbott Tower 
Dayton, OH 45402 

Max Young 

Director of Advanced Engineering 
Utilimaster Corp. 

BOX 585 
Wakarusa , IN 


46573 
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COMPOSITE MATERIALS AND STRUCTURES PROGRAM 

rpi workshop on Composite Materials 
and Structures for Rotorcraft 


Rensselaer Polytechnic Institute 
Rotorcraft Technology center 

September 10 A 11, 1986 


Attendees 


Mr. M. Adam 

Messerscmitt - Bolkow - Blohm Gmbh 
Poetfach 801160 
Munich 80 800 

Federal Republic of Germany 


Prof. Marco Borri 

Cntr for Advancement of Computational Mechanics 
Georgia Institute of Technology 
School of Civil Engr 
Atlanta, ga 30332 


Mr. G. Reis Alsailler 

S601 Oakmont Lane 

Fort Worth, TX 76112 


Dr . Douglas Cairns 

Massachusetts Institute of Technology 

Cambridge, Ma.' 02139 


Dr . Gary Anderson 
US Army Research Office 
P0 Box 12211 

Research Triangle Park, NC 27709 

Mr. Robert Arden 
3580 Royal Hills Drive 
St. Louis, M0 63129 

Mr . Boas Avitzur 

watervliet Arsenal 

watervliet, NY 12189-4050 

Mr. D. J. Baker 

us Army Aerostructures Directorate 
NASA Langley Research center 
Hampton, VA 23665 

Dr. OUvler Bauchau 
Dept, of Mechanical Engineering 
Rensselaer Polytechnic institute 
Troy, NY 12181 


Dr. Raymond G. Carlson 
Project Engineer - RiD 
Sikorsky Aircraft 
N. Main. St. 

Stratford , CT 06602 

Dr. Robert Cheatham 
Sikorsky Aircraft 
N. Main St. 

Stratford , CT 06602 

Mr. Milton Chin 

Area Representative 

AMOCO Performance Products Inc. 

1300 Lakeside Ave 

Cleveland, OH 44114 

Dr. Inderjit Chopra 
University of Maryland 
College Park, MD 20742 

Dr. Tsu-Wei Chou 
University of Delaware 
Newark, Delaware 19711 
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Dr. Mark Darlow 
Dept, of Mechanical Engineering 
Rensselaer Polytechnic Institute 
Troy, NY 12181 

Dr . M . Degener 
Institute of Aeroelasticity 
DFVLR 

Bunsenstrasse 10 
D - 3400 Goettingen 
Federal Republic of Germany 

Mr. H. Benson Dexter 
Sr. Researcher 

NASA Langley Research center 
Hampton, VA 23665 

Dr. R. Judd Diefendorf 
Department of Materials Engineering 
Rensselaer Polytechnic Institute 
Troy, New York 12181 

Mr . Michael Dubberley 
Attn: AIR 5302 
Dept, of Navy 
NAIR Systems Command 
Washington, DC 20332 

Dr. John Dugundji 

Massachusetts Institute of Technology 
Cambridge, Ma. 02139 

Dr. Wolf Elber 

Chief, Materials and Structures Area 
US Army AATD 
Hampton, VA 23665 

Mr. Gary Farley 

US Army Aerostructures Directorate 
NASA Langley Research Center 
Hampton, VA 23665 

Mr. J. C. Fish 
University of Maryland 
College Park, MD 20742 

Mr . Michael Frengley 
MDHC 

4645 South Ash Avenue 
Tempe, Az. 85282 


Dr. Joseph Goldberg 
MS S529A 
Sikorsky Aircraft 
N. Main St. 

Stratford, CT 06602 
Dr . Dan Good 

Aviation Applied Technology Directorate 
SAURT-TY-ATS 

Ft. EustiS, VA 23604-5577 

Mr. Kenneth Grina 
Vice President - Engineering 
Boeing vertol Company 
PO Box 16858 

Philadelphia, PA 19142 

Maj . G . Haritoe 
Program Manager 
AF0SR/NA 
Bolling AFB 

Washington, DC 20332 

Mr. Derek S. Hart 
Chief, structures Design 
Boeing Vertol Company 
PO Box 16858 
Philadelphia, -PA 19142 

Mr. Irwin J. Kenigsberg 
MS S330A 
Sikorsky Aircraft 
N. Main St. 

Stratford , CT 06602 

Dr. Paul A. Lagace 
Rm 33-313 

Massachusetts Institute of Technology 
77 Mass Ave 

Cambridge, MA 02139 
Ms . Renee Lake 

US Army Aerostructures Directorate 
NASA Langley Research Center 
MS 340 

Hampton, VA 23665 

Prof. S. w. Lee 
University of Maryland 
College Park, MD 20742 

Mr. L. L. Levine 
Structure Technology 
Sikorsky Aircraft 
N. Main St. 

Stratford, CT 06602 


Mr. Samuel P. Garbo 
45 High Hill circle 
Madison , CT 06443 



Dr. Robert G. Loewy 
Institute Professor 
Rensselaer Polyechnic Institute 
Troy, New York 12181 

Mr. Giancarlo Maffioli 
Agusta SPA 
Via Luini IS 
Aibizzate 210U7, Italy 

Mr. Azar P. Majidi 
University of Delaware 
Newark, Delaware 19711 

Mr. Leonard Marchinski 
6 East Avenue 

Mount Carmel, PA 17851 

Dr . Roger Massey 
Vice Pres. • Engineering 
Kaman Aerospace Corp. 

Blue Hills Ave, POB 2 
Bloomfield, CT 06002 

Mr. Bill Matthews 

Army Materials Technology Laboratory 
Watertown, MA. 

Mr. James A. Hondo 
Project Manager - Helo 
E. I. duPont Nemours Co, Inc 
Centre Road Bldg 
Wilmington, DE 19898 

Mr . Robert Murrill 
Sr. Project Engr. 

McDonnell Douglas Helicopter Co. 
4645 South Ash Avenue 
Tempe, AR 85282 

Mr. Mark Nixon 

US Army Ae restructures Directorate 
NASA Langley Research Center 
MS 340 

Hampton, VA 23665 
Mr. Peter C. Ogle 

Mgr, Composite A/P Design a Struct. 
Sikorsky Aircraft 
N. Main St. 

Stratford, CT 06602 


Dr. Kevin O'Brien 
US Army Ae restructures Directorate 
NASA Langley Research Center 
MS 188E ARTA 

Hampton, VA 23665-5225 

Mr. David J. Parker 
Sect. Spvr - Comp. MAP Engr 
Sikorsky Aircraft 
N. Main St. 

Stratford , CT 06602 
Mr . K . Pfeifer 

Messerscmitt - Bolkow - Blohm Gmbh 
Poetfach 801160 
Munich 80 800 

Federal Republic of Germany 

Mr . Robert Pinckney 
Boeing Vertol Company 
P0 Box 168S8 

Philadelphia, PA 19142 

Mr. Pietro Putelli 
Agusta SPA 
Via Luini 15 
Aibizzate 210U7 , Italy 

Dr. Lawrence w. Rehfield 
Aerospace Engineering 
Georgia Institute of Technology 
Atlanta, ga 30332 

Dr . Dale Reisdorfer 
Bell Helicopter TEXTRON 
PO Box 482 

Port Worth, TX 76101 

Mr. John W. Roman 
Deputy Manager, Adv. Composites 
Grumman Aircraft 
MS B10-2S 

Bethpage, LI, NY 11714 

Mr. Glenn Rossi 
MS P32-38 

Boeing Vertol Company 
PO Box 168S8 

Philadelphia, PA 19142 


Mr. Satesh salpekar 
c/o Dr. Kevin O'Brien 
NASA tangley Research center 
MS 188 E ARTA 

Hampton, VA 23665-5225 

Dr. A. A. Salzberg 
University of Maryland 
college Part, MD 20 742 

Dr. A. D. Stemple 
University of Maryland 
College Part, MD 20742 

Dr . Iradj Tad jbakhsh 
Department of Civil Engineering 
Rensselaer Polytechnic Institute 
Troy, New York 12181 

Dr. T. Vizzini 
University of Maryland 
College Part, MD 20 742 

Dr . Steven Wlnckler 
Dept, of Mechanical Engineering 
Rensselaer Polytechnic Institute 
Troy, New York 12181 
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PERSONNEL 


C o -Principal Investigators 

Loewy , Robert G. , Ph.D. 

Wiberley , Stephen E. , Ph.D. 

Senior Investigators 

Bauchau , o . , Ph.D. 

(Structural dynamics , advanced 
composites )* 

Diefendorf * , R. J. , Ph.D. 
(Fabrication, resin matrix, fiber 
behavior, interfaces)* 

Dvorak*, G. , Ph.D. 

(Metal matrix composites, damage, 
fracture & fatigue of composites)* 

Krempl, E. , Dr.Ing. 

(Fatigue studies, failure criteria)* 


Shephard, M. S., Ph.D. 

(Computer graphics, finite element 
methods )* 


Sternstein*, S. S., Ph.D. 

(Failure analysis, matrix behavior, 
moisture effects)* 

Research Staff 


Institute Professor 

Professor and Chairman of Chemistry 


Assistant Professor of 
Aeronautical Engineering 


Professor of Materials 
Engineering 


Professor and Chairman of Civil 
Engineering 


Professor of Mechanical Engineering, 
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8:15-8:45 
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8:45-10:15 


10:15-10:30 
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10:30-12:30 


m WORKSHOP ON COMPOSITE MATERIALS 
AW STRUCTURES FOR ROTORCRAFT 

FINAL PROGRAM SCHEDULE 
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Aircraft, Stratford, Connecticut. 

‘"Application of Composites to the Multi Service V- 22 Osprey", 
V. Hart, Boeing Vertol , Philadelphia, Pennsylvania. 
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"Crash Energy Absorbing Composite Helicopter Structure", 

G. L. F anley, Aerostructures Directorate, AVSCOM, Langley 
Research Center, Hampton, Virginia. 

3:45-4:00 BREAK 

Session IV : CHAIRMAN - R obert Pinckney ( Boeing Vertol) 

Generic Structural Elements 


4:00-6:00 


6:00-6:30 

6:30 


"Evaluation of Composite Components on the Bell 206L 
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